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The electronic spectrum of a two-dimensional square lattice in a perpendicular magnetic field has 
become known as the Hofstadter butterfly [Hofstadter, Phys. Rev. B 14, 2239 (1976)]. We have 
calculated quasi-one-dimensional analogs of the Hofstadter butterfly for carbon nanotubes (CNTs). 
For the case of single- wall CNTs, it is straightforward to implement magnetic fields parallel to the 
tube axis by means of zone folding in the graphene reciprocal lattice. We have also studied perpen- 
dicular magnetic fields which, in contrast to the parallel case, lead to a much richer, pseudofractal 
spectrum. Moreover, we have investigated magnetic fields piercing double-wall CNTs and found 
strong signatures of interwall interaction in the resulting Hofstadter butterfly spectrum, which can 
be understood with the help of a minimal model. Ubiquitous to all perpendicular magnetic field 
spectra is the presence cusp catastrophes at specific values of energy and magnetic field. Resolving 
the density of states along the tube circumference allows recognition of the snake states already 
predicted for nonuniform magnetic fields in the two-dimensional electron gas. An analytic model of 
the magnetic spectrum of electrons on a cylindrical surface is used to explain some of the results. 

PACS numbers: 73.63.Fg,73.22.-f,73.43.-f,73.43.Qt 



I. INTRODUCTION 

The availability of new materials for nanoelectronic 
research allows for a detailed test of the emergence of 
the quantum physical nature of electrons, via transport 
or optical measurements. Carbon nanotubes^'^'^'^'^'^ 
(CNTs) are an example of a very peculiar electronic 
material, due to the extreme confinement of electrons 
on their 7r-conjugated "walls". In these systems, many 
mesoscopic phenomena such as single-electron charging,^ 
and conductance quantization,^ as well as effects typical 
for semiconductor physics like s-like excitons,^^ can be 
observed already at room temperature. 

Since the prediction of band structure effects of car- 
bon nanotubes in parallel external fields by Ajiki and 
Ando in 1993,^^ it took only a few years until clear hall- 
marks of a single quantum flux being tethered within a 
tube section were found experimentally in optical^^ and 
transport^^'^^ measurements. For magnetic fields per- 
pendicular to the CNT axis, theoretical predictions were 
made shortly after, first using a perturbative approach 
around the Fermi energy, and later also using a tight- 
binding model^SiiS Only recently, a first experimen- 
tally accessible effect of perpendicular magnetic fields — 
anomalous magnetoconductance — was predicted^^ and 
observed^ A very similar effect for strong electric fields 
has also been found by numerical studies^^ and has yet to 
be confirmed experimentally. The use of magnetic fields 
to further investigate the interplay between elastic mean 
free path, phase coherent length, and electron-electron 
interaction was also successfully adopted. ^^'^^'^^ 

From the purely theoretical perspective, carbon nan- 
otubes in strong perpendicular magnetic fields repre- 
sent a very interesting case of study. Closely related to 
graphene, their energy spectrum shows strong similar- 
ities with that of the two-dimensional (2D) honeycomb 



lattice which again forms a variation of the fractal 
butterfiylike pattern discovered by Hofstadter^ in 1976 
and studied intensely since that time from various points 



of 



27.28.29.30.31.32 



Yet the quasi- ID nature and the 



curvature of CNTs set their energy spectra clearly apart 
from the fractal and perfectly periodic images obtained 
in 2D lattices. 

In this paper, we will describe a method of comput- 
ing and visualizing the spectrum of carbon nanotubes 
(for a prototypical example see Fig. P). This method 
will be demonstrated on a number of single- and double- 
wall CNTs (SWCNTs and DWCNTs) of different chiral- 
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FIG. 1: (Color online) Density of states of a (6,6) CNT in 
dependence on an external magnetic field parallel (left) or 
perpendicular (right) to the tube axis. For every value of the 
magnetic field, the DOS is unity normalized over energy. The 
units = $o/r^7r and Bq = ^oMpiaquette (see text) are 
scaled such that the physical field scale is the same for both 
segments of the plot. 
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FIG. 2: (Color online) The structure of a CNT [here, a (3,3)- 
CNT]: the hexagonal lattice of a graphene sheet is rolled up 
in such a way that the chiral vector (n, m) becomes the cir- 
cumference of the resulting cylinder. Magnetic fields paral- 
lel to the tube axis pierce the tube cross section r^vr, while 
perpendicular magnetic fields pierce the wall made up from 
hexagonal plaquettes. 

ity and diameter. The study of the local distribution of 
the spectral density will shed some light on the relation 
between the spectrum of a planar sheet of graphene and 
that of a CNT, strongly affected by curvature and finite 
size. A closer look at the spectrum will reveal the pres- 
ence of cusp catastrophes, which are closely related to 
the quenching of the Bloch state velocity, induced by a 
magnetic field. 

For magnetic fields parallel to the tube axis, the nat- 
ural unit is that of one flux quantum per tube cross 
section r^yr (see Fig. |2j. For a general (n, m) CNT 
the tube radius can be obtained with simple geomet- 
rical arguments 27rr = V3m^ + 3n^ + Smndcc^ where 
dec = 1-42 A denotes the carbon-carbon distance. This 
immediately gives the parallel magnetic field Bq needed 
to pierce one flux quantum $o = h/e through an (n^m) 
CNT. For perpendicular magnetic fields, the scale is ruled 
by the field necessary to enclose a flux through a sin- 
gle benzene ring, the plaquette of grap hene and carbon 
nanotubes of area Apiaquette = 3^3/4 die ^ 5.24 A^. 
Because of this extremely small area we obtain Bq = 
*^*o /^plaquette = 79 X 10^ T, which is, of course, out of 
experimental reach. It is straightforward to get the re- 
lation between the parallel and perpendicular field scales 
as 



For a typical SWCNT with ~1 nm diameter, this gives 
a value of Bq ^ 5 x 10^ T. It is thus understandable 
that multiwall CNTs (MWCNTs) present a more inter- 
esting object for magnetic field experiments: For a typ- 
ical MWCNT with a diameter of 20 nm, as a matter of 
fact, one can already observe the first Aharonov-Bohm 
oscillations accessible at around 12 T parallel fields. 



As shown in this work, however, even for perpendicular 
fields low-field signatures could be visible within experi- 
mentally accessible field ranges if one takes into account 
the external shell of a MWCNT. 

This paper is organized as follows. We first give defi- 
nitions, introduce the method of computation and visu- 
alization, and point toward general features observable 
in quasi- ID systems. In Sec. Ill, we then do a system- 
atic study of SWCNTs, including an analytic model and 
a detailed view of the range of experimentally accessi- 
ble fields. In Sec. IV, we proceed with an analysis of 
the effects of the interwall interaction in DWCNTs on 
the magnetic spectrum and introduce a minimal model, 
closing with a discussion of the results in the last secion. 

II. DEFINITIONS, METHODS, AND 
OBSERVABLES 

Lattice electrons in arbitrary external magnetic 
fields. All numerical calculations in this work are based 
on a tight binding Hamiltonian of the form 

where the indices denote the atomic orbitals. For the 
single-orbital approximation used hereafter, these coin- 
cide with the label of the atom so that H can be repre- 
sented by means of the matrix elements Idn.rj between 
TT orbitals centered on the atom i and j at the position 
Vi (vj) of the CNT molecular network. 

An external magnetic field is implemented using the 
Peierls substitution>24 Based on the principle of mini- 
mal coupling p p — eA, the effect of a magnetic 
field B = rot A is absorbed in the translation opera- 
tor T (R) = exp \{p — eA) • R. In the tight-binding 
Hamiltonian, this is reflected by representing the hop- 
ping matrix elements jij = (^i|7Y|^j) between two tt 
orbitals and localized at sites Vi and rj as 

7,, (B) = 7° exp (^i^ £ dr ■ As (r)) . (2) 

The bare hopping at zero magnetic field acquires a 
complex phase expressed as an integral along the bond 
direction dij = rj — Vi . 

With the CNTs oriented parallel to the z axis, it is 
advantageous to choose a gauge in such a way that Ab 
is independent of r^. This is provided, e.g., by 

AB{r) = (0,r,5||,r^5x), (3) 

giving a magnetic field B = (5x,0,5||) with known 
components perpendicular and parallel to the tube axis. 
Throughout this work fields will be consider either per- 
pendicular {B\\ = 0) or parallel {B^ = 0) to the tube axis. 
Arbitrary angles are of course possible as well, showing 
the expected crossover of both regimes. 
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FIG. 3: (Color online) Scheme to illustrate the physical meaning of the butterfly plots. An external magnetic field distorts 
the band structure of a CNT in an intricate way. For any fixed magnetic field, the DOS and the transmission can be obtained 
directly from the band structure, (a), (b), and (c) are sections of the two right panels of the DOS and transmission vs E and 



Having chosen a linear gauge further simplifies the in- 
tegration in Eq. to a product: 



lij exp 



•2^^ A 
i—dij • Ab 



^0 



(4) 



In the presence of a perpendicular magnetic field, it is 
thus necessary to consider the exact coordinates of the 
molecular structure at hand rather than — as sufficient for 
parallel or vanishing magnetic fields — their simple topo- 
logical connectivity. 

Density of states. For such an r;^-independent 
gauge field, the Hamiltonian of any quasi- ID periodic 
structure like a CNT stays periodic in the presence of a 
magnetic field. This allows the use of the Bloch theorem 
to derive the corresponding band structure. As can be 
seen in Fig. O the band structure is in general strongly 
distorted by an applied magnetic field. The density of 
states (DOS) can be determined from the magnetic band 
structure Ei^ik^B) via 



IS 



the band index and 



where 

3dccVm^ + + mn/ gcd (3n, n — m) the length of 
the unit cell of an (n^m) CNT. Ei){k^B) is obtained 
by direct diagonalization of the CNT Hamiltonian 
via the Bloch ansatz (see Appendix Since we 

work in a basis of one orbital per atom, the number 
of bands TVb equals the number of atoms in the unit 



cell N 



4(r 



nm) I gcd (3n, n — m)^ growing 



with the diameter and dependent on the helicity angle 
(deviation from the armchair or zigzag configuration). 
The plot of the DOS directly reflects the distortions of 
the band structure caused by the magnetic field. 



Alternatively, Green-function-based approaches allow 
one to resolve the DOS within different atoms in the same 
unit cell by introducing the local density of states 

PLDOS.. {E, B) = --Im a.,,., {E, B) , (6) 

TT 

where Gn.n is the space-diagonal component of the lat- 
tice Green function matrix (see Appendix IB|) 



g {E, B)=[E- n{B) + iO+] ' , 



(7) 



Of course by tracing the LDOS within the different atoms 
of the same unit cell, one can restore the full DOS as 
Pbos{E, B) = jj: Xlr, Pldos^^ {E, B). 

Butterfly plots. To capture the continuous evolu- 
tion of the band structure with growing magnetic fields, 
it is very convenient to visualize the DOS in butterfly 
plots, as illustrated in Fig. |31 The resemblance to the 
well-known Hofstadter butterfly of 2D lattice electrons^) 
becomes very clear for CNTs of large diameter (see Fig.[7| 
below) . A common feature to butterfly plots of all quasi- 
ID systems are the pronounced band edges, caused by 
van Hove singularities in the DOS.^^ 

In Fig.^ a (6,6) CNT Hofstadter butterfly is plotted as 
a reference for further comparisons. For the parallel field, 
the behavior is perfectly periodic for integer multiples of 
the flux quantum — h/e penetrating the tube cross 
section r^7r. Starting as a metallic CNT at 5 = 0, the 
gap opens and closes periodically.^^ 

For perpendicular fields with their natural scale of one 
flux quantum per graphene plaquette, the overall behav- 
ior is not periodic. This can be understood due to the 
presence of plaquettes at various angles toward the field, 
capturing different, in general incommensurate, fractions 
of the flux quantum. However, a number of features from 
the underlying graphene structure are still visible at the 
diameter-independent scale of Bq . 
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FIG. 4: (Color online) DOS in a graphene ribbons of infinite length and various widths and internal orientations, pierced 
perpendicularly by magnetic fields. Each ribbon can be classified as an unrolled CNT: The "chiral" vectors refer to the 
SWCNT which, when unrolled, would result in the corresponding planar ribbon. The density of states is normalized to the 
number of atoms per unit cell to give a comparable visual appearance. 



Important to note is the difference in the behavior for 
small fields: while the parallel field causes a linear Zee- 
man split of the states with opposite angular momentum, 
small perpendicular fields generally cause quadratic en- 
ergy shifts. 

All plots are of course symmetric in the magnetic field 
sign, which is why only half butterfly plots are shown. 
Different is the case of the E ^ —E symmetry which 
is related to the particle-hole symmetry. The latter is 
present in the 7r-orbital description of SWCNTs but is 
broken by the interwall interaction in DWCNTs. 

Transport observables in quasi- ID systems. As 

can also be seen in Fig. 01 it is straightforward to apply 
the same scheme not only to the density of states, but just 
as well to other properties like the quantum mechanical 
transmission T of a quasi- ID system. The latter is the di- 
mensionless zero-temperature conductance after the Lan- 
dauer theory of phase-coherent transport:^^ G = GkT^ 
where Gk = 2e^//i is the conductance quantum and in- 
verse of the von Klitzing resistance. The calculation of 



the transmission, which involves a renormalization pro- 
cedure for the semi-infinite carbon nanotube leads^^ by 
means of the energy-dependent injection rates Fl/r^ and 
the Green function Q projected on a finite nanotube par- 
tition, can be cast into the Fisher and Lee formula^ 



Still, for a periodic structure — as is the case for the 
systems at hand — the quantum mechanical transmission 
is simply a band-counting algorithm, and as such con- 
tains less information than the band structure itself or 
the DOS. This is very different from magnetotransport 
through finite CNTs: Scattering at the contacts leads 
to resonant tunneling, resulting in spectroscopy of the 
electronic states of the finite tube.^^'^^ This spectrum 
may show strong dependence on magnetic fields, even 
in regions of flat bands, resulting in quantum-dot-like 
physics^ 
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Relation to 2D periodic structures. It is impor- 
tant to note some similarities, but also some fundamen- 
tal differences between the butterfly plots of quasi- ID 
structures and those in the original work by Hofstadtepi 
and later generalizations^^-^^ which handled 2D periodic 
structures. Starting out from an analogous Hamiltonian 
and also using the Peierls substitution to implement the 
magnetic field, the most striking difference is that, for 
a 2D periodic structure, it is not possible in general to 
choose a gauge in such a way that the resulting Hamilto- 
nian has the same translational symmetry as the under- 
lying system. For rational values of the magnetic flux per 
unit cell, one can still find a larger effective unit cell but, 
for irrational values, this is not possible at all, which 
ultimately leads to the fractal structure of the energy 
spectrum found by Hofstadter, similar to that displayed 
in the lower panel of Fig. [7| In contrast, the quasi-lD 
structure of CNTs results in a fixed number of bands, 
leading to a pseudofractal spectrum, with the recursion 
of self-similarity limited by the transverse length-cutoff 
of the system. 



Graphene ribbons. Since the recent experimental 
success in isolating single sheets of graphene^Si^ the ex- 
otic Dirac-like electronic structure has become the focus 
of several studies. Epitaxially grown graphene has been 
used to laterally confine electrons and determine coher- 
ence lengths studying weak-localization effects in magne- 
totransport measurements.^^ For understanding the re- 
lation between the butterfly of a 2D graphene sheet and 
these quasi- ID carbon nanotubes, it is instructive to take 
a look at graphene ribbons as an intermediate step. An 
(n, m) graphene ribbon is simply a planar "unrolled" 
(n, m) CNT, periodic in one dimension and finite in the 
other. As in the original Hofstadter butterfly,^ the ribbon 
butterfly plots are periodic as a function of the perpen- 
dicular magnetic field due to the equal flux piercing any 
hexagonal plaquette forming the honeycomb lattice (see 
Fig.^J. As the ribbon width increases the butterfly plots 
tend to the Hofstadter butterfly of a graphene layer as 
visible in the bottom panel of Fig. [7| 



Cusp catastrophes. One striking detail ubiquitous 
in butterfly plots are the cusp catastrophes appearing 
at specific positions of energy and magnetic fields (see 
Fig. 0). These are points where, with changing mag- 
netic field, some band is continuously deformed from a 
strictly monotonic curve into a band with two adjacent 
zero-group- velocity points. At the exact point where this 
mathematical catastrophe happens, both the first and 
the second derivatives of E{k) are zero. A wave packet 
of this energy and momentum will have both its velocity 
and its spreading suppressed leading to a special kind of 
localization not unlike that of Landau levels. 




A (Si) 

FIG. 5: (Color online) Cusp catastrophes are ubiquitous in 
butterfly plots. The band structure at magnetic fields below, 
at, and above the critical magnetic field shows the smooth 
transition from a strictly monotone band into a third-order 
parabolic band with changing magnetic field. 



III. SINGLE- WALL CARBON NANOTUBES 

For SWCNTs we consider only the radial p orbitals — 
forming the most electronically relevant tt bands — and 
only interactions between nearest neighbors, setting 
7ij = 7o = 2.66 eV. This has been shown to be an ex- 
cellent approximation in explaining electronic structural 
and transport properties of SWCNTs. ^'^'^ The on-site 
energy Si = Eq is constant for all atoms and defines the 
Fermi energy = sq of di neutral CNT. Ignoring an 
offset in the energy, we can simply choose Sq = 0. Zee- 
man splitting could also easily be included in this cal- 
culation as So = ±^/iB^/2 and would result visually in 
an overlay of two butterfly plots sheared against each 
other linearly with growing magnetic fields. The inten- 
sity of this effect at the critical plaquette field scale is 

^/^B ^plaquette = 9.1 cV. 

The special case of parallel magnetic fields: 
Shortcut via the zone-folding method. As an al- 
ternative to calculating the electronic bands of a SWCNT 
via the procedure described above, one could calculate 
the spectrum of graphene and then apply periodic bound- 
ary conditions in the angular direction of the CNT {zone 
folding). For magnetic fields parallel to the tube axis, 
this method is still applicable: the phase gathered by an 
electron moving on a closed loop around the tube cir- 
cumference can simply be included in the boundary con- 
ditions. This results in a shift of the allowed discretized 
quasimomenta in the reciprocal space. For perpendicular 
magnetic fields, however, this method breaks down and 
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FIG. 6: (Color online) Two semiconducting SWCNTs of sim- 
ilar diameter as the (6,6) CNT in Fig. Q] The band gap os- 
cillates irregularly with increasing perpendicular field. The 
large unit cell of the the chiral (6,5) tube leads to bands with 
low dispersion as soon as the rotational symmetry is broken 
by the perpendicular magnetic field. 



one has to consider the full geometry of the CNT. 



A. Structural properties 

Chirality dependence. Several features can be 
found when comparing the magnetic spectra of tubes 
with different chiralities though similar diameter (see 
Figs, n and (i) The behavior of the gap around 
the charge neutrality point E = is very helicity 
dependent: a parallel magnetic field always opens and 
closes the gap periodically as a consequence of the in- 
teger number of fluxes per nanotube cross section. This 
phenomenon is independent of whether the tube is metal- 
lic or semiconducting at 5 = 0. In contrast, for per- 
pendicular fields there are distinctions: Armchair CNTs 
stay strictly metallic for any perpendicular field, as can 
be understood from supersymmetry arguments^ On the 
other hand, the gaps of the two semiconducting CNTs in 
Fig. El do open and close in an aperiodic, though oscilla- 
tory, pattern. The gap closes to zero in single points of 
specific values of Bj_ and opens again. Closer observa- 



tions of a larger set of CNTs reveal that this also happens 
for semimetallic tubes like the (3n, 0) zigzag CNTs. The 
precise opening and closing pattern carries an intrinsic 
complexity; its statistical behavior, however, seems to 
depend mostly on the number of atoms in the unit cell, 
(ii) Another general effect of the large unit cell in the 
(6, 5) CNT is that in this chiral tube with its large num- 
ber of plaquettes at different angles towards the field, 
the high symmetry of the original system is broken down 
very efficiently by the magnetic field, resulting in a larger 
number of bands of very low dispersion. The magnetic 
field effectively localizes the electrons in nonpropagating 
Landau-like states. 

Diameter dependence. Figure [7| shows the evo- 
lution toward the graphene Hofstadter butterfly of the 
magnetic spectrum of armchair CNTs as a function of 
their diameter. The (200,200) SWCNT has a diameter 
of 27 nm, comparable to the external shell of a typi- 
cal MWCNT, and is thus of great interest. The visual 
effect — resembling watercolors "flowing" toward the right 
hand side — can also be understood by a simple picture. 
For large enough diameters, the CNT consists of regions 
of nearly flat graphene, each at a different angle toward 
the magnetic field, thereby experiencing a different nor- 
mal component of the magnetic field, as visible in the in- 
set of the middle panel of Fig.[7| Since the DOS is an av- 
erage over the LDOSs at the different unit cell atoms, one 
ends up with a sum of different graphene Hofstadter but- 
terflies, stretched to different effective fields, the stretch 
being minimal where the magnetic field is normal to the 
tube wall and maximal where it is tangential. Overlaying 
these differently stretched graphene butterflies results in 
the "flowing" appearance of the butterfly of large diam- 
eter tubes. 

At the lower and upper energy edges, one can clearly 
see the emergence of linear Landau levels and the char- 
acteristic fractal structure of the graphene butterfly is 
unmistakenly visible at the same scale of the magnetic 
field. In fact near the top and bottom of the graphene 
TT energy band (of width 2W = 67 = 16 eV), electrons 
have an effective mass of m* = 2Ji? /S^d^Q ^ 0.95me, 
leading to a cyclotron frequency uj{B±) = eB^/m*^ so 
one could write 

E = TW±huj{B)(^n^^^ (8) 

with n = 0,l,2,..., which fits nicely with the numerical 
results (as indicated by the straight lines in the bottom 
panel of Fig.^J. 

Around the Fermi energy, the Dirac-like dispersion 
of graphene leads to the so-called relativistic Landau 
levels, ^^'^^ following 

E = ±VY^2nheB±_ (9) 

with n = 0,1,2,... and the Fermi velocity = 
Z'^dcc/2h. These levels, which can be clearly observed 
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FIG. 7: (Color online) Comparing different diameters: the 
large (200,200) CNT bears strong resemblance to the Hof- 
stadter butterfly of graphene (Ref. 25«) combined with the 
curvature effects (details in text). The straight lines at the 
lower left corner of the graphene butterfly (bottom panel) in- 
dicate the Landau states obtained from an effective mass con- 
tinuum theory [see Eq. The parabolic lines near in 
the same plot indicate the relativistic Landau levels obtained 
from the Dirac-like dispersion of graphene [see Eq. @]. 



in the Hofstadter butterfly of graphene (left edge in 
the bottom panel of Fig. [TJ, are also responsible for 
the recently observed anomalous quantum Hall effect of 
graphene 424^ 

Snake states. The view of the total DOS of a large 
SWCNT as the sum of different contributions from the 
regions at various angles around the tube circumference 
can be confirmed by taking a look at the LDOS at in- 
dividual atom positions. As can be seen in Fig. [HI the 
LDOS at ^ = 0, where the magnetic field pierces the wall 
perpendicularly, resembles very much the butterfly of the 
planar graphene sheet. The electrons here show very low 
dispersion, similar to Landau levels. At = 7r/2, on 
the other hand, the magnetic field is tangential to the 
CNT wall and therefore has far less effect on the electron 
dispersion. An understanding of the electronic states in 
these regions can be gained by considering classical elec- 
trons confined to the surface of a cylinder: As the effec- 
tive magnetic field (the projection of the field onto the 
tube normal) changes sign at = 7r/2, the curvature of 
an electron trajectory will also switch orientation each 
time the electron crosses this "equator" line, leading to 
a snakelike movement of the electron^Sii^ 



B. Analytical model 

In order to shed some light of intuition on our results, 
we may consider the physics of a structureless hollow 
cylinder, a tubule, in a perpendicular magnetic field. 
(Similar systems in parallel magnetic fields have been 
studied before42iS) This system bears some similarity 
to a "Hall bar", with the crucial difference that it does 
not have borders that could carry edge states. Instead, 
it has two flanks where the magnetic field is tangential 
to the tube and therefore the radial component of the 
magnetic field — which is the effective field experienced 
by electrons confined to the cylinder surface — vanishes. 
To understand where charges do accumulate, we consider 
this system in cylindrical coordinates (6>, z) at fixed ra- 
dius r. By Eq. a perpendicular magnetic field leads 
to a gauge field 

A{0) = B^r sine 

in cylindrical coordinates. With this, the Hamiltonian of 
an electron restricted to the tube surface becomes 



1 



2mr^ 



pI 



eB±r sini 



which can be viewed as that of an electron in 2D with 
periodic boundaries in a nonuniform magnetic field. 
A similar system — a 2D strip ranging over [— L/2,L/2] 
in the y direction and infinite in the z direction, placed 
in a linearly varying magnetic field B = BoyCx — was 
first studied in 1992 by Miiller,^^ who identified two new 
classes of states: one at finite magnetic field propagating 
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FIG. 8: (Color online) Decomposition of the density of states into the contributions of particular atoms (identified by their 
angle $ toward the magnetic field direction). A plaquette at angle captures a flux of B^Apiaquette cos^. The region at $ = 
experiences a perpendicular field piercing the tube wall, very much as in the plain graphene sheet (Fig. The regions at 
$ = 7r/2 experience a field tangential to the tube wall, leading to a much smaller flux per plaquette, resulting a the stretched 
impression of the butterfly. The DOS butterfly over the whole CNT unit cell is an overlay of these and many intermediate 
pictures. In the angle-resolved plot for B = 0.2 Bq one can see a smooth transition between a region with Landau levels and 
a region with normal band dispersion. For the stronger field, the systems goes through two oscillations along the angle. 



perpendicularly to the field gradient direction with loop- 
ing trajectory and low velocity, the other around the line 
5 = 0, propagating in the opposite direction at higher 
velocity with a snakelike trajectory. To solve our system, 
we can exploit the z invariance and do an ansatz for the 
wave function: ^ (^, z) = ipk^ iO) e^^^^. Our problem re- 
duces to that of a particle in one dimension with a kz 
dependent potential: 



Vkz (0) 



2mr^ 

VkA^) = —{hkz-eB^rsmef. (10) 
2m 

For \hkz\ < \eB±r\^ this potential has two minima at 
^min = 7r/2 ± 8ii ccos {Hkz / eB ±r) . A harmonic approxi- 
mation at either of these minima yields the approximate 



Hamiltonian 



2 

Pe 



1 

with the spectrum 

n 



[eB^rf - [hk^f 



Enih) = — J {eB^rf -{hk, f{n+ 1/2). 
mr ^ 



From this dispersion relation, we can directly retrieve the 
group velocity 



(n + 1/2) hkz 



'{eB±r) 



(11) 



{hkz 



The wave functions in the harmonic potential are lo- 
cated around the minima ^min, so for low energies we 
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can say in reverse that at each angle we find pre- 
dominantly electrons with the longitudinal wave vector 
kz{0) = {eB±r/h)sinO. Placing this into Eq. (|TT|) . we 
can retrieve an expression for the velocity of electrons 
moving at certain angles: 

Vn(0) = — ^ (n + 1/2) tan l9. 
mr 

Now, the divergency at 6 = ± (7r/2) originates from the 
fact that the harmonic approximation breaks down when 
the two minima of Vk^ {0) meet at this angle. Apart from 
this, however, one can see clearly the angular separation 
of electrons moving in both directions and the localiza- 
tion in Landau-like states at ^ = and ^ = tt, where the 
magnetic field pierces the tube wall normally (see Fig.[3|). 

Such a continuum model can only be expected to hold 
for the CNTs at low magnetic fields with 5±Apiaquette 

As it turns out, some of the effects visible at 
higher fields can be understood qualitatively by study- 
ing a model of intermediate complexity: a square lat- 
tice cylindrical tube of lattice constant a. Coming from 
the continuum model and following Ref. [H^, we can 
replace the continuous coordinates by integer indices: 
{am^an) := {rO^z). Using a tight-binding model with 
on-site energy Eq and hopping parameter 70, the Hamil- 
tonian acts on a wave function in the following way: 

(m— l,n) ~r ^(m+l,n)/ 

- 70 (e-*^(™) + e"^("*) *(™,„+i)) 

where the phase factor (m) = (27rar5^/^o) sin (ma/r) 
originates from the Peierls substitution Eq. As in the 
continuum, the invariance in the z direction can be ex- 
ploited, now using a Bloch ansatz due to the discreteness 
of the system: 

^(m,n) - e %)m • 

This leads to a finite Hamiltonian for any fixed kz G 
[— 7r/a, 7r/a): 

Hkzi^m = -70 (V^m-l + V^m+l) + Vk, (m) l/jm, 

xr / N fi raB± ma\ 

Vk, (m) = £0 - 270 cos kza - 27r— — sin . 

V r J 

The most significant difference to the effective potential 
of the continuum model Eq. is the replacement of 
the square law by a cosine one. This has the effect that 
the potential does not grow indefinitely for large mag- 
netic fields, but instead oscillates, forming several min- 
ima at various angles ^, as seen in Fig. El In combination 
with the discretization of the angle, this potential leads 
to the formation of a complex pattern in the angular de- 
pendence of the density of states, as it can be observed is 
the LDOS at high magnetic fields also displayed in Fig. [SI 




FIG. 9: (Color online) Analytical solution of the continuum 
model. Top panel: the A;^ -dependent effective potential of free 
electrons confined to a continuum cylinder in a magnetic field 
perpendicular to the tube. Middle panel: first three eigen- 
states of the harmonic approximation to the above potential 
for various kz (highlighted wave functions correspond to the 
selected potentials in the first panel). Each wave function 
is shifted to the corresponding energy. Superimposed are the 
lines followed by the extrema of the wave functions. The same 
pattern can be found in the top left panel of Fig.|Hl where the 
maxima of the DOS show the maxima of the various energy 
eigenstates. Bottom panel: The /c^-dependent effective poten- 
tial of a discretized tube showing a large number of minima. 
States located in narrow minima have higher energy, so the 
low-energy spectrum is mainly determined by the widest po- 
tential minima. 

To capture more details in a model, an appropriate 
step would be the implementation of the correct disper- 
sion at the Fermi energy: The characteristic cones at 
the Fermi points of graphene can be approximated by a 
Dirac-like Hamiltonian. For a detailed study of the mag- 
netic spectrum of Dirac-electrons on a cylindrical surface, 
see ^45^54] . 

C. Experimentally accessible perpendicular 
magnetic fields 

In recent experiments, optical transitions in CNTs 
were studied in magnetic fields up to 75 T.^^ New ex- 
periments are in preparation to go up to 200 T and even 
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FIG. 10: (Color online) Zoom into the butterfly of two differ- 
ent armchair SWCNTs. The scales, including the color scale, 
are chosen according to the scaling law given in the text to 
produce comparable representations of the data. The tubes 
correspond to diameters of 2.7 nm (top) and 27 nm (bottom). 
The inset in the bottom panel illustrates the shape of the peak 
at the Fermi level. The white lines crossing the plot in the 
upper panel are caused by small avoided crossings in the band 
structure. 



2-3 kTiS 

In perpendicular fields of this magnitude, as displayed 
in Fig. 03 the first onset of the band structure distortions 
can be seen clearly in large CNTs, comparable with the 
outer shell of typical MWCNTs, measuring up to tens of 
nm in diameter. 

In particular the plots show very clearly the rapidly 
changing van Hove singularities, resembling those of 
Fig.nfoi" a (6,6) CNT, but at much lower magnetic field 
scale. Moreover, the E = graphene state, which is 
due to the peculiar distortion of the Dirac-like linear dis- 
persion into a strongly nonlinear one,^^ emerges at lower 
fields with increasing diameters. 

Most notable is the scaling law that can be found in 
the butterfly plot of large tubes at low fields near the 
Fermi energy: For two different armchair CNTs with the 



chiral vectors (m, m) and (m^, m') it can be expressed as 



This scaling is followed approximately already for small 
CNTs and becomes very precise for large diameters, con- 
verging toward a DOS that is reproducible from a model 
of Dirac electrons on a continuum cylinder. The peak at 
the Fermi energy also follows this scaling law. Within the 
region of scaling, the maximum of the peak at £^ = Ey 
grows exponentially with the magnetic field while its in- 
tegral grows linearly. 

It is important to note that the scaling is not an effect 
of the curvature, but of the discretization of the transver- 
sal momentum, since it can be observed in graphene rib- 
bons as well. 



IV. DOUBLE- WALL CARBON NANOTUBES 

While SWCNTs and MWCNTs have been studied in- 
tensely over the past 15 years, it has only recently be- 
come possible to produce DWCNTs of high purity and 
qualityr^iS fueling the interest in details about the in- 
terwall interaction. Previous studies have shown an in- 
teresting interplay between magnetic fields parallel to 
the DWCNT axis and the interwall interaction near the 
Fermi energy. A minimal Hamiltonian of a DWCNT 
can be set up as 



n 



E 



by defining the intrawall interaction as described for 
SWCNTs. For the interwall interaction, we can fix the 
hopping coefficients as 



{B) = f3 cos i^ij exp 



X exp 



1— d., • Ab 



3.34 A, 6 = 0.45 A, and i^ij and 



.^J stand for the angle and the absolute distance between 



where p = 70/8, a 

the two TT orbitals {{i^j)) centered at positions Vi and r 
belonging to two different shells. ^^'^^ 

As a representative example, the butterfly of a 
(6,6)@(11,11) DWCNT is displayed in Fig. [ill For the 
parallel magnetic field, the interwall interaction gives rise 
to features at two levels: looking at individual bands, one 
can observe van Hove singularities crossing and some- 
times avoid a crossing. The complex rules whether a 
crossing is avoided are not obvious from studying the 
butterfly plot only, but can be deduced by looking at 
the band structure and taking into account the various 
symmetries of the system. At larger scale in the butter- 
fly plot, one finds a modulation of the pattern crossing 
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^ i(6,6)(a)(ll,ll) DWCNT without intershell interaction 




|(6,6)@(11,11) DWCNT with intersliell interaction | 
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FIG. 11: (Color online) Butterfly plot of a (6,6)@(11,11) 
double- wall CNT. In the upper panel, the interwall interac- 
tion is switched off, resulting in an overlay of the butterflies of 
two independent SWCNTs. In the lower panel, the interwall 
interaction gives rise to a number of new features (see text 
for details). 



from E = - SeV at B\\ = to E = E^ ^ SeV at 

To understand this phenomenon, we have studied a 
single unit cell in a magnetic field perpendicular to the 
plane of the resulting concentric ring (Fig. [T^ . The 
spectrum shows a periodic behavior of the modulation 
with a period ^^^terwaii ^ $o/(^^ccMnterwaii/4) relat- 
ing to the area of the minimal circular path between 
both walls. This period and the shape of the modula- 
tion are independent of the diameter of the DWCNT. 
A physical explanation for this modulation is as follows. 
Like the interaction in a two-atomic molecule, splitting 
two atomic orbitals into a bonding and an antibonding 
molecular orbital, the interwall interaction may also hy- 
bridize SWCNT states of the same energy into bonding 
and antibonding DWCNT states. However, the system 
has an approximate rotational symmetry, so the interwall 
interaction may only hybridize states of the same angular 
momentum. 

At zero magnetic field, the angular momentum of the 
states at the bottom of the spectrum is zero in both shells. 
This allows hybridization, causing a split in the hybrid 
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FIG. 12: (Color online) Upper panel: Spectrum of a system 
of two concentric atomic rings. Atom spacing and coupling 
inside each ring are taken from graphene. The distance 5r 
between the rings as well as the the parametrization of the 
coupling between the rings follow those given in the text for 
DWCNTs. The sketch displays the prevalent links between 
the shells. Even though the geometry is irregular, the area 
of circular paths is very near to integer multiples of dccSr/2, 
leading to a clear periodicity of the modulation in the spec- 
trum. Lower panel: An isolated unit cell of a DWCNT with 
the same radii as the planar double ring of the upper panel. 
This system has smallest closed loops at an angle against 
the magnetic field, resulting in an effective smallest area of 
dccSr/4: and leading to a doubling of the period. Further- 
more, the system has two atoms in the rotational periodic 
cell, leading to two interlaced modulations. 



spectrum. At the upper end of the spectrum, the angular 
momentum does not match, prohibiting a hybridization. 
This is the cause for the strong electron-hole asymmetry 
visible in DWCNT butterfly plots. 

By switching on a parallel magnetic field, the effec- 
tive angular momentum is shifted by the Aharonov-Bohm 
phase gathered on a circular path around the tube. This 
shift depends on the cross section of the path, so it is dif- 
ferent for the two shells. Therefore, the energy at which 
both angular momenta match depends on the magnetic 
field, causing the splitting region to travel over the en- 
ergy range, which leads to the visible modulation in the 
parallel field butterfly plots. 

For a simplified model — a double-wall square lattice 
tube — the modulation does follow a single cosine-shaped 
curve, as shown in the upper panel of Fig. El In compar- 
ison, the DWCNT shows an additional complexity: the 
underlying honeycomb lattice of graphene has a unit cell 
containing two atoms, resulting in two intertwined cosine 
curves, the second just becoming visible at the edge of 
Fig.ini 

For fields perpendicular to the axis of a DWCNT, the 
only large-scale effect caused by the interwall interac- 
tion observable in the butterfly plot is the hybridization- 
induced splitting already described for zero field. With 
growing field, this effect disappears, and the plot shows 
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no remarkable global patterns. 



V. CONCLUSIONS 

The magnetic spectrum of two-dimensional infinite lat- 
tice electrons gives rise to the well-known Hofstadter but- 
terfly. In this paper, we have shown that quasi-one- 
dimensional lattice electrons exhibit a spectrum which 
does resemble the fractal structure of the Hofstadter but- 
terfly but with a flnite cutoff due to the transversal con- 
flnement. We have calculated such pseudofractals for 
carbon nanotubes, a material at the focus of many nano- 
electronic studies also in relation to the presence of ex- 
ternal magnetic flelds. We have calculated the density of 
states (butterfly plots) of several single wall carbon nan- 
otubes and we could show (i) the strong dependence of 
the magnetic spectrum on the underlying chiral indices; 
(ii) the emergence of the graphene Hofstadter butterfly 
at increasing nanotube diameter. In particular, perpen- 
dicular flelds induce an aperiodic and pseudofractal mag- 
netic spectrum. Periodic structures have been obtained 
for graphene ribbons, demonstrating that the aperiod- 
icity of the perpendicular fleld butterfly plots is due to 
the incommensurability of the magnetic flux captured by 
elementary (hexagon) plaquettes of a CNT oriented at 
different angles towards the external fleld. 

By studying the angle-resolved electronic structure of 
a SWCNT one can observe the emergence of snake states 
already predicted for nonuniform magnetic flelds in a Hall 
bar.^^ In our case, we have been able to devise an ana- 
lytical model for the states at the top and the bottom 
of the energy spectrum by means of an effective mass 
approximation. In this latter case a continuum theory 
can capture the striping of the wave function along the 
region of the tube with zero normal fleld. Inversely, near 
the Fermi level, one cannot bypass the Dirac neutrino 
nature of the electronic states. We have interpreted the 
wave function striping by writing a Harper equation^^ for 
square lattice electrons with a cylindrical geometry. 

While the effects of parallel flelds are of comparably 
simple nature in SWCNTs (being an expression of the 
Aharonov-Bohm oscillations due to a rigid shift of the 
graphene band structure sampled via the zone-folding 
method), this is not the case for DWCNTs. The electron- 
hole symmetry of tt bands in SWCNTs is broken once two 
shells are put in interaction. The resulting hybridization 
of inner and outer states could be clearly understood by 
means of two interacting Aharonov-Bohm rings. 

Experimentally relevant effects have been calculated 
for SWCNTs of diameter of typical external shells in 
MWCNTs. There underlying multifractal structure like 
that of Hofstadter can be observed already at a few tens 
of tesla, and an outstanding scaling law for the DOS at 
low magnetic flelds near the Fermi energy has been given. 
The latter applies also to graphene ribbons and is intrin- 
sically related to the massless dispersion at the charge 
neutrality point. 



This study, though systematic, could not include very 
interesting issues which also deserve careful investigation, 
such as the effects of disorder on the butterfly plots of 
SWCNTs. Disordered SWCNTs can be thought in fact 
a model for the external shell of MWCNTs. More atom- 
istically one could study the influence of the interwall 
interaction of the structure of large diameter DWCNTs 
(also reasonable models for MWCNTs) 3 
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APPENDIX A: HISTOGRAM METHOD 

A histogram method is the simplest method to get the 
magnetic spectrum of a quasi-lD system. It is also very 
efficient if the complete energy range has to be calculated. 
Starting from a periodic Hamiltonian of the form: 



H 



V 



Hi Ho Hi 

Hi Ho Hi 



I 



(Al) 



one can use Bloch theorem to get an effective Hamilto- 
nian: 



'H, 



^H\ 



where a is the length of the unit cell. Numerically scan- 
ning the ID Brillouin zone — 7r/a < k ^ i\ j a with a 
uniform distribution, one can now diagonalize the flnite 
matrix i^eff (/^) for each value k. The resulting eigenval- 
ues from this diagonalization are counted in a linear his- 
togram over the full energy range and normalized to the 
total number of states. Depending on the resolution of 
the k sampling, this histogram will become an arbitrarily 
good approximation to the density of states. Figures 
ini[ZI(CNT panels), 21 and ITTl were calculated using this 
method. 

The calculation of the data in Fig. [TUIwas heavily op- 
timized by using an adaptive k sampling in combination 
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with a linear interpolation to reduce the number of di- 
agonalizations in regions of smooth band structure and 
increase the precision at band edges. 

APPENDIX B: GREEN FUNCTION METHOD 

Another, more flexible method is that using Green 
functions: The bulk Green function Q (E) of the infinite 
CNT can be calculated very efficiently by the following 
method. 

The periodic Hamiltonian in Eq. (|A1|) is used as start- 
ing point of a recursive decimation scheme: 

hI°\e) = Ho, 

Hi°\E) = Hi 

With each recursion, the length of the effective unit 
cell is now doubled by decimating out every second cell: 

^01 - ^01 7' '^01 ^ 

rr(n+l)/^x _ {n) (n)TT{n) 

^10 y^) - ^10 7' '^10 ' 

where 7^^^ = {E ^vq — Hq^^)~^ and 77 is a small positive 
numerical value, chosen smaller than the desired energy 
resolution but large enough to provide fast convergence 
and numerical stability. 

Convergence is reached for n > if \\H^'1\e)\\ + 

Ili^Q^ ^ (^) II < ^ foi" some matrix norm || • || and some small 
cutoff e. We can then retrieve the bulk Green function 
from the converged Hq as: 



With the original Hamiltonian (|A1|) expressed in a tt- 
orbital tight-binding basis, the resulting Green function 
Sbuik is a matrix in the same atomic basis of one unit 
cell. Therefore, the local density of states in each atom 
is directly given by 



PLDOS.(^) = -^Im[a(^)J, 



summing up to the poos {E) = J^li Pldos^ {E). In the 
same run, the surface Green functions Q^^^ (E) can be 
used to calculate the transmission through the system 
using the Fisher-Lee relation^^ with a single unit cell se- 
lected as conductor, as shown in Fig. 

The Green function method and the histogram method 
both give numerical approximations to the same mathe- 
matical quantity, but their numerical errors are very dif- 
ferent: while the former method tends to give ffuctua- 
tions that show up as grainy structure in ffat areas of the 
butterffy plot, the latter suffers from sampling problems 
around van Hove singularities. Both errors have to be 
countered with very high resolution scanning and down 
sampling of the data. The data presented in the figures of 
this article typically took several hours to weeks of com- 
putation time on standard PCs [Intel(R) Pentium(R) 4, 
3 GHz]. 

The work presented here was done using the fol- 
lowing Open Source (R)^^ software: Python as pro- 
gramming language, NumPy (Refs. 'g? and and 
SciPy (Ref. 69) for numerical computations, PyTables 
for data storage and handlingfffi matplotlib for data 
visualization,^^ inkscape for figure preparation,^^ and 
TeXmacs for authoring. 
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